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1 Introduction 



A superintegrable system is a system that is integrable (in the sense of Liouville- Arnold) 
and that, in addition to this, possesses more constants of motion than degrees of freedom. 
At this point we must note that the maximum number N of functionally independent 
constants of motion for a system in a ci-dimensional manifold is N = d — 1. Thus if a 
Lagrangian (or Hamiltonian) system has n degrees of freedom then, as the phase space is 
2n- dimensional, we have that the maximum number of independent constants of motion 
is N = 2n — 1. There are three well known examples of this very particular class of sys- 
tems, namely, the Kepler problem, the isotropic harmonic oscillator, and the non-isotropic 
oscillator with commensurable frequencies. The two-dimensional harmonic oscillator is 
a system trivially integrable since it can be considered as a kind of "direct sum" of two 
systems with one degree of freedom. If the oscillator is isotropic then it has the angular 
momentum as an additional integral of motion. If the oscillator is non-isotropic the an- 
gular momentum is not preserved as the potential is not central; nevertheless when the 
quotient of the two frequencies is a rational number then the system has another addi- 
tional integral. Concerning the three-dimensional Kepler problem, it possesses not only 
the energy and the angular momentum as constants of motion, but also the Runge-Lenz 
vector; only five of these seven integrals are functionally independent since in this case 
the number of degrees of freedom is n = 3. In these three cases it is well known that all 
the orbits became closed for the case of bounded motions. This high degree of regularity 
(existence of periodic motions) is a consequence of the superintegrable character. 

Fris et al [1] studied the two-dimensional Euclidean systems which admit separability in 
two different coordinate systems and obtained four families of potentials V r , r = a,b,c,d, 
possessing three functionally independent integrals of the motion. A very important 
point is that these three constants of motion were linear or quadratic in the velocities 
(momenta). In fact, if we call superseparable a system that admits Hamilton- Jacobi 
separation of variables (Schrodinger in the quantum case) in more than one coordinate 
system, then quadratic superintegrability (i.e., superintegrability with linear or quadratic 
constants of motion) can be considered as a property arising from superseparability. The 
first two families 

v a = y (x 2 + y 2 ) + - 2 + - 2 , 

V b = |w 2 (4x 2 + y 2 ) + k 2 x + -§ , 

can be considered as the more general deformations (with strengths k 2 , k$) of the 1:1 and 
2:1 harmonic oscillators preserving quadratic superintegrability (the other two families, 
V c and Vd, were related with the Kepler problem). The superintegrability of V a was later 
on studied by Evans [21 [3] in the more general case of n degrees of freedom. 
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A natural generalization of V a is given by the following potential 

V a (ui, w 2 ) = \(u\3? + w 2 V) + % + 

x y 

that contains a more general harmonic oscillator with anisotropy. This new potential is 
only separable in Cartesian coordinates and its superintegrability (if it exists) must be of 
higher-order ('higher-order superintegrability' means that some of the integrals of motion 
are polynomials in the momenta of order higher than 2). Therefore, the method of the 
multiple separability cannot be used and it must be studied by making use of a different 
approach. In fact, the superintegrability of this nonlinear system, for the case of rational 
ratio of the frequencies, was first proved by Ranada et al in [4] using the properties of 
the isotonic (or singular) oscillator and the Pinney-Ermakov equation [5J for obtaining 
a complex factorization. More recently Evans et al [6] and Rodriguez et al (7[ E] have 
also proved this property using the geometric formalism of dimensional reduction. It have 
been proved that certain nonlinear integrable systems may arise as reductions of very 
simple systems defined in higher- dimensional spaces (see e.g. [9j [10] ). In this particular 
case the authors start with an harmonic oscillator in a higher-dimensional space and then 
the dimensional reduction introduces the nonlinearities but in a way that preserves the 
superintegrability. We note that this geometric method has been also used for proving 
the superintegrability of the Kepler-Coulomb system with nonlinear terms [8j [Til [12] . 

Now, we present a new method to prove the higher-order superintegrability of this 
nonlinear system. This new method, that is more straightforward than the previously 
known methods, is directly related with the approach presented in [I] (but without making 
use of the properties of the Pinney-Ermakov equation) and, at the same time, it is also 
related with some of the results obtained in JTJ [8] . 

The plan of the article is as follows: in Sec. 2 we present the problem from the La- 
grangian viewpoint and we study the superintegrability of the two dimensional potential 
V a (ui, U2) using the existence of a complex factorization as an approach. First we consider 
the linear harmonic oscillator and then we prove the superintegrability of the nonlinear 
system. The constants of motion of the nonlinear system are explicitly obtained. In 
Sec. 3 we analyze the general case of n degrees of freedom. It is written by using the 
Hamiltonian approach and it has a more generic character. We start with a very general 
Hamiltonian H and then we introduce the restrictions to obtain superintegrability (but 
without making use of the property of separability). The idea is that this more general 
approach could be used in future papers as a starting point for the search of other more 
general superintegrable systems. Finally in Sec. 4 we make some comments and present 
some open questions. 
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2 Complex factorization and superintegrability 



2.1 Superintegrability of the linear Harmonic Oscillator 

The two-dimensional harmonic oscillator 

T 1 / 2 , 2\ 1 / 2 2, 2 2\ 

Lho = 2 ( v x + v y) ~ 2(^1 x +u 2 y ) 

has the two partial (one-degree of freedom) energies, I± = E x and I 2 = E y , as fundamental 
integrals. The superintegrability of the rational case, oj\ = n x u , u 2 = n y u , with integers 
ti x , fiy , can be proved by making use of a complex formalism [T51 [EH [15] • Let K h i — x,y, 
be the following two complex functions 

K x = v x + i n x u x , K y = v y + i n y u y , 

then we have the following time-evolution 

d d 

—K x = i n x u K x , —K y = i n y u K y . 

Thus, the functions defined as 

K,, . i,j = x,y, 

are constants of motion. The two real functions l-fGxl 2 and \K yy \ 2 are proportional to the 
energies E x and E y and concerning K xy , since it is a complex function, it determines not 
one but two real first integrals, lm(K xy ) and Re(K xy ). So, we have obtained four integrals 
but, since the system is two-dimensional, only three of them can be independent. We can 
choose I\ = E x , I 2 = E y , and ^3 = lm(K xy ) as the set of fundamental constants of motion 
(the other constant I± = R,e(K xy ) can be expressed as a function of E x , E y , and Im^^)). 

As an example, for the Isotropic case, u\ = u 2 = co , we obtain 

h = — lm{K xy ) = xv y - yv x , 
h = Re(K xy ) = v x v y + io 2 xy , 

(^3 is the angular momentum and I4 the component F xy of the Fradkin tensor [16]) and 
for the first non- isotropic case, U\ = 2uq, u> 2 = u>o, we arrive to 

I3 = — lm(K xy ) = (xv y - yv x )v y - u 2 xy 2 , 
h = Re(K xy ) = v x v 2 y + u 2 (Axv y - yv x ) y . 

The integral -Z3 of the 3:1 oscillator will be cubic and, in the general the 
function J 3 will be a polynomial in the velocities (momenta) of degree n x + n y — 1 . 
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2.2 Superintegrability of the nonlinear system 



The technique presented in the previous section proves not only the super-integrability 
of the rational case but also the existence of a complex factorization for the additional 
constant of motion. Now we assume, as starting point of our approach, that if a new 
system can be obtained by a deformation of the harmonic oscillator then it must be also 
endowed with a similar property. 

The analysis will proceed in two steps. 

Step 1. The associated quadratic equation 

Let K = a + ib be such that dK/dt = inuK. Then the function K 2 defined as 
K 2 = K 2 = (a 2 — b 2 ) + 2\ab satisfies dK 2 /dt = 2inuK 2 . Thus the time evolution of the 
functions 

K 2x = {v 2 x - n 2 x u 2 x 2 ) + 2in x u xv x , K 2y = {v 2 - n 2 y uly 2 ) + 2\n y u yv y , 
is given by 

— K 2x = 2 i n x u K 2x , — K 2y = 2 i n y u K 2y , 
and hence the complex functions K 2i j defined as 

K 2ij = (K 2i p (X 2 *.y\ i,j = x,y, 

are constants of motion. 

The two complex functions, Kij and K 2i j, must be considered as two alternative ways 
to prove superintegrability but, of course, the first one is simpler than the quadratic. As 
an example we have \K 2xx \ = \K XX \ 2 and \K 2yy \ = \K yy \ 2 . In the general 
the functions /3 = lm(K 2xy ) and I4 = Ke(K 2xy ) will be a polynomials in the velocities 
(momenta) of degree 2(n x + n y — 1) — 1 and 2(n x + n y ). So, if the study is restricted 
to the harmonic oscillator, it is better to use the function since it leads to simpler 
expressions for the constants of the motion. 

Step 2. Introducing a deformation in the "quadratic equation" 

Let us now consider the following (F, G)-dependent family of potentials 

V(n x , ny, F, G) = \u 2 (n 2 x x 2 + n\y 2 ) + \F{x) + \G{y) . 

The problem is to determine which particular values of the functions F and G can preserve 
the existence of a complex factorization. At this point we assume that, in order to solve 
this problem, it is more convenient to use the quadratic equation since it seems as more 
'deformable' than the linear one. 
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Let us denote by Aj and Bj, j = x,y, the following functions 

A x = v 2 x - n 2 x u 2 x 2 + F(x) , B x = 2n x u xv x , 
A y = v 2 y - n 2 y uly 2 + G{y) , B y = 2n y u yv y . 

Then, if we require for the functions Aj and Bj a time-evolution of the form 

d n d 

-A 3 = -2n 3 u Q B 3 , - 



Aj = - 2njU Bj , —Bj = 2n j w Aj , j = x,y, 



we arrive (we omit the details) to the first-order differential equations 

xF' + 2F = , yG' + 2G = , 

with solutions F = k\/x 2 and G = k 2 /y 2 with arbitrary constants hi and k 2 . Therefore, 
only if F and G have this particular structure, the complex functions M 3 - defined as 
Mj = Aj + i Bj play, in this nonlinear case, a similar role to the complex functions Kj of 
the linear case. 

Proposition 1 Consider the non-linear potential 

V a (n x , n y ) = \u 2 {n 2 x x 2 + n 2 y y 2 ) + ^ + ^ 

representing an harmonic oscillator with rational ratio of frequencies, oj\ = n x uo, u 2 = 
n y u , and inversely quadratic nonlinearities and let us denote by Mj, j = x,y, the follow- 
ing two complex functions 



M x = (vl - n 2 x u 2 x 2 + ^ + 2in x u xv x , M y = (v 2 - n 2 y uly 2 + ^ + 2in y u yv y . 

Then, the complex functions M i3 defined as 

M l3 = (M t p (M;p , i,j = x,y, 
are constants of the motion. 

Firstly, the moduli of M x and M y are given by 

\M X \ 2 = A(E 2 X - hnlu 2 ) , \M y \ 2 = A(E 2 - k y n 2 y uj 2 ) 
The time-evolution of the functions M x and M y is given by 

4- M x = 2 i n^wo M x , ^- M y = 2 i n„w M„ . 
at at 
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Thus we have 

j t M xy = n y {M x f n ^ (M* y ) n *M x + n x (M x ) ny (M^"*" 1 ) M * 
= {2iu ){n y n x - n x n y )(M x ) ny (M*) Ux = . 

As in the linear case, M xy can be considered as coupling the two degrees of freedom. 

Hence the potential V a (n x ,n y ) is superintegrable for any rational value of the the quo- 
tient U2/0J1. In the particular 1:1 case, u\ = co 2 = Uq, the potential reduces to the V a 
potential 

y a (l,l) = y a = !„ V + ! / 2 ) + ! 5 + ^ 

(that is the only superseparable potential in this family) and in the general case, u>\ = 
n x (jj , u 2 = n y u , the potential V a (n x ,n y ) represents a generalized V a potential with a 
non-isotropic n x :n y oscillator (note that in the 2:1 case, the potential V a (2, 1) is different 
from the family V&). 

In the general the functions Ke(M xy ) and lm(M xy ) will be polynomials in 

the velocities (momenta) of degree 2(n x + n y ) and 2(n x + n y ) — 1 respectively. The real 
part takes the form 

Re(M xy ) = 2^ +n y\E x ) n y(E y ) n * + Xu 2 J 3 , 

where A is a numerical coefficient. Thus the additional constant of V a (n x ,n y ) is in fact 
the function J3 that is of degree 2(n x + n y — l). If we denote by -Z3 the constant of motion 
of the associated linear system (of degree n x + n y — 1), then the additional third integral 
J3 can be written as follows 

j 3 = n + h4 w) + k 2 jf l) + ki4 20) + k x k 2 j^ + ... + k^4 0nx) ■ 

That is, the integral J3 of the non-linear system (hi 7^ 0, k 2 7^ 0) appears as a deformation, 
not of the function J3 itself, but of its square If (this property was already mentioned in 

ID)- 

Next, we give the expressions of the constant J3 for the three first cases: 

(i) Potential V a (l, 1) corresponding to a central (isotropic) harmonic oscillator. 

In this case we have u\ = u 2 = ojq and the integral of motion J 3 of the associated 
linear system is just the angular momentum, J 3 = (l/uo)lm(K X y) = xv y — yv x . 
Then we have 

Re(M xy ) = AE x E y - 2w 2 J 3 , 

with J3 given by 
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(ii) Potential V a (2,l) corresponding to a non-isotropic 2:1 oscillator with frequencies 
coi = 2co , co 2 = wo- 
lf we denote by J 3 the integral of motion of the associated linear system 

h = 77— lm ( K xy) = {xv y - yv x )v y - u 2 xy 2 , 

then we obtain 

Re(M xy ) = 8E X E 2 - 8u 2 J 3 , 

with J 3 given by 



h = J 3 + i^) v l + ^(y^ - 2xv v) 2 + 



2y^ yux ^ Vy> ^ 2x* 

(iii) Potential V a (3, 1) corresponding to a non-isotropic 3:1 oscillator with frequencies 
coi = 3co , oo 2 = wo- 
lf we denote by J3 the constant of motion of the associated linear system 

h = — lm(K xy ) = 3{xv y - yv x )Vy + ul(yv x - 9xv y )y 2 , 

then we obtain 

Re(M xy ) = 16E x E 3 y - 2u 2 J 3 , 

with J 3 given by 

1/3 — 1% -\- k\ 1/3 ^ -|- J3 ^ k\k2 1/3 ' ^ "I - 1/3 ' ^ + k±k 2 t/3 ' ^ -|- /c^ t/3 ' ^ j 
with the functions J 3 10 \ <4 01 \ -T^ 1 '^, and Jg 0,3 ' 1 , given by 

10 2/ 2 01 3 

T (i,i) 12-u y T (o,2) 3 2 

j(i,2) = _3_ j(0 ,3) = 9x 2 ^ 

3 a; 2 ?/ 2 ' 3 

Summarizing, within the V a (n x ,n y ) family, only in the particular isotropic 1:1 case 
the function J 3 is quadratic in the velocities and, because of this, only in this case the 
superintegrability arises from separability in two different coordinate systems. In all the 
remaining cases, V a (n x ,n y ) is a superintegrable but not superseparable potential. 
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3 Hamiltonian formalism and n degrees of freedom 



The previous section was directly focused on the two dimensional potential V a (n x ,n y ). 
Now we present the study of the general case of n degrees of freedom but using a rather 
different approach. The idea is to start with an integrable but very general Hamiltonian H 
and then look for the properties to be satisfied by H in order to admit superintegrability 
determined by appropriate complex functions K^. 

Let H be the following Hamiltonian 

n n 

i=i i=i 

defined in a 2n-dimensional phase space T*Q (Q is the n-dimensional configuration space) 
endowed with the standard Poisson bracket 



_ /dR dS _dR dS\ 
1 ^ \dxi dpi dpi dxi' 



It is clear that H is integrable; now we set out the problem to determine the expressions of 
the functions Vi{ Xi ) to admit superintegrability (but without making use of the property 
of separability). 

Let us now define a set of n functions (possibly complex) linear in the momenta 

Ki = Pi + fi(xi) , i = l,2,...,ra, 
(note that we have chosen the coefficient of pi equal to 1). Then, it is evident that 

{K l ,K j } = 0, 

but 

{K i ,H} = fl(x i ) Pi -V!(x i )^0, 

(except in the trivial case fi(xi) = constant and Vi(xj) = constant). Our next step is to 
write the Hamiltonian in terms of the Ki functions. Since we assume that these functions 
can be complex, we impose 



i=i 



where K* is the complex conjugate of Ki. So we have 

i^i 2 = p-+p,(/(^) + r(^)) + i/(^)i 2 . 
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If we assume that the Hamiltonian H is quadratic in the momenta, and without linear 
terms, then the functions should be pure imaginary. So we arrive to 

n 

Ki = pi + i gi(xi) , V = \^2 / g i {x i ) 2 , 

i=i 

with gi(xi) real functions. 

3.1 From integrability to superintegrability via the functions K{ 

The Hamiltonian H is trivially integrable, the n functions \Ki\ 2 are functionally indepen- 
dent constants of motion and the Hamiltonian is half their sum. In order to study its 
superintegrability, we make use of the tensor Ty defined by 

n 

T l3 = K t K] with H= 1 -J2Tn. 

i=i 

Since the Poisson bracket of the functions Ki with the Hamiltonian H is now given by 

{Ki,H} =g' i (xi)(ip i -g i (x i )) = ig'^x^Ki, 

then the off diagonal tensor components have the following Poisson bracket with the 
Hamiltonian 

{% H} = {K h H}K* + Ki{K*, H} = i (g^) - g'^xj))^ . 

Thus if the g'^Xj) are a numerical constant (independent of i) then the functions Tij are 
constants of the motion. This happens with 

gi(xi) = u Xi , i = l,...,n, 

so we get the isotropic harmonic oscillator in n dimensions which is certainly superinte- 
grable: 

n n 
i=l i=l 

3.2 Generalizing the Hamiltonian Hi 

If we consider, instead of the tensor T, a new tensor T with components 

f l3 =K^(K*r>, i,j = l,...,n, 
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(where n« are positive integers) we can repeat the arguments above and obtain 

{%#} = n j {K i ,H}K?-\K*)^ + n i K?{K*,H}{K*r- 1 
= i(n j g^(x i )-n i g' j (x j ))T i:i . 

Thus if the functions <7j(xj) satisfy the relations 

n j9i( x i) ~ n i9j( x j) = , i,j = l,...,n, 

then the functions (components of the new tensor T) are constants of the motion. 
In this way we obtain the Hamiltonian of the rational anisotropic harmonic oscillator 
(non-central harmonic oscillator with rational ratio of frequencies) 

n n 
ZJ 1 \ ^ 2,1 2 \ ^ 2 2 
H ^ = 2 l^Pi + 2 U }^ n i X i 

i=l i=l 

which represents the most simple superintegrable generalization of the Hamiltonian H 1 . 



3.3 Generalizing the Hamiltonian H2 

In order to obtain a superintegrable generalization of the Hamiltonian H 2 , we first consider 
the squares of the functions K i 



K t = Pi ~ n \u^x\ + 2 i riiUoXiPi , i = 1 , . . . n, 



satisfying 



{Kl H 2 } = {K h H 2 }Ki + Ki{Ki, H 2 } = 2i mu> K? , 
and then we introduce a deformation of the real part 

M i = Pi - n \uQx\ + hi(xi) + 2iniU XiPi , i = 1,. . .n. 

Now let us consider the following Hamiltonian 

n n n 



i=i 



i=i 



i=i 



where Vi(xj) are some functions to be determined by imposing that |Mj| 2 Poisson com- 
mutes with H. We have 



{\Mi\ 2 , H} = 2(h'i - 2VD p\ + 2 |2(/iJ - 2^% - n,Wa\ ? (4/ii + + 2^F/) 
so that we arrive at 

h'i - 2V( = , 4/i, + + 2xiV! = , 



Pi, 
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with solution 

1 k- 
Vi(xi) = -hi{xi) = i = l,...,n, 

(up to inessential additive constants). Hence, the above Hamiltonian H becomes 

n n n , 

jj 1 \^ 2,1 2 \ 22,1 \ ki 

i=l i=l t=l * 

in such a way that the n functions |Mj| 2 coincide with the square of the n energies E$ up 
to a constant 

|M| 2 = 4(E 2 - hnW) , Ei = \p 2 + \^ 2 n 2 x 2 + A . 

The important point is that the functions Mj satisfy 

{M i ,i2g} = 2ira i w M < . 
Hence, if we denote by the functions defined by the products 

/ k ■ \ n j 

M,j = Ml h (M*) n > = U - n 2 tu 2 x 2 + 4 + 2 i mu; x iPi ) 

2 2 2 2 , o • \ 

- // ; ^o ./' ; • yj ~ ZirijUoXjPjj , 

then we have 

{M ij ,H 3 } = 0, i,j = l,...,n, 

what means that both the real part and the imaginary part of My are constants of the 
motion for H 3 (the diagonal functions Ma are real and, as we have seen, they correspond 
to the energies Ei). As the system has n degrees of freedom (and the phase space is 2n- 
dimensional) the maximum number of independent constants of motion is N = 2n — 1; so 
it is clear that not all of these quantities will be functionally independent but, nevertheless, 
we can extract from them a fundamental set of 2n — 1 functionally independent invariants. 
So we conclude that the Hamiltonian H 3 is superintegrable for all the values of the integer 
numbers n» and for arbitrary values of the constants ki. 

Next we prove the above statement (it is always possible to extract, from the large 
number of constants of motion, a fundamental set of 2n — 1 functionally independent 
functions) by following the same arguments used for the harmonic oscillator (isotropic 
and nonisotropic). We first recall the following two points 

(i) In the isotropic n-dimensional case, the Fradkin tensor F [16] is represented by a 
symmetric n-dimensional matrix F^ = piPj + uo 2 XiXj, i,j = 1,2, ... ,n, so that 
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it provides a total set of (l/2)n(n + 1) constants of motion. The integrability is 
consequence of the n diagonal entries F# (related with the energies Ei) that are 
independent in a trivial way; so we have 



For proving the superintegrability we recall that every nondiagonal function F^ only 
depends of the four variables (xi,Pi, Xj,pj). Thus we can add, for example, the n—1 
entries i^j+i of the upper-next-diagonal so that we obtain 

dF n AdF 22 A . . . AdF nn AdF 12 AdF 23 A . . . AdF n „ ln ^ . 

Hence, these N = 2n — 1 constants of motion are functionally independent. 

(ii) In the non-isotropic n-dimensional case, the method of the complex factorization [131 
E] discussed in subsection (12. ID leads to a complex Hermitian n-dimensional matrix 
Kiji hj — 1)2, • • - ,n, with Kij only depending of the four variables (xi,Pi,Xj,pj); 
so the property of independence 

dK u AdK 22 A . . . AdK nn Ad(lmK 12 )Ad(lmK 23 )A . . . Ad(lmK n ^ ln ) ^ , 

is also true in this case (of course in this complex case it is also possible to choose 
the Real part of the functions K^). 

Now, in the case of the nonlinear Hamiltonian H 3 , we have obtained a complex Her- 
mitian n- dimensional matrix My, i,j = 1,2, ... ,n, that can be considered as a nonlin- 
ear deformation of the matrix K^. But, since only depends of the four variables 
(xi,Pi,Xj,Pj), we also have the following property 

dM u AdM 22 A . . . AdM nn Ad(lmM 12 )Ad(lmM 23 )A . . . Ad(lmM n _ ln ) ^ . 

Thus, we have proved the existence of a set of iV = 2n — 1 functionally independent 
constants of motion. 

We close this section with an interesting property. If we return (for ease of notation) 
to two degrees of freedom, then in the case of the harmonic oscillator {k\ — 0, k 2 = 0), 
the constant of the motion I4 = Ke(K xy ) can be obtained (up to a factor) as the Poisson 
bracket of I3 = lm(K xy ) with the energy E x (for example, in the 1:1 case the component 
F xy of the Fradkin tensor arises as the Poisson bracket of the angular momentum J with 
E x ). This property is preserved by the deformation (hi 7^ 0, k 2 ^ 0) and remains true for 
the nonlinear system. 

Proposition 2 The Poisson brackets ofKe(M xy ) andlm(M xy ) with E x are given by 



dFuAdF 22 A . . . AdF nn + . 



{lm(M xy ) , E x } 
{Re(M xy ) , E x } 



2u} n x n y Re(M xy ) , 
- 2u n x n y Im(M xy ) . 
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Note that the second Poisson bracket means that lm(M xy ) (degree 2(n x + n y ) — 1) is 
just the Poisson bracket of J3 with E x 

{J 3 ,E x } = -^Llm(M xy ) 

where A is a numerical coefficient. 

4 Comments and open questions 

As stated in the introduction the superintegrability of V a (n x ,n y ) was firstly proved in jl] 
and then in [5J [7J [5] by the use of different methods. Of course these methods are all 
correct (the dimensional reduction has been previously applied to the study of a certain 
number of integrable systems) but we think that the approach presented in this paper 
(deformation of the quadratic equation) has the great advantage of possessing a great 
level of elegance and simplicity. Moreover it is related with one of the more fundamental 
properties of the harmonic oscillator. 

We note that although this method is rather different from the the dimensional reduction 
they have in common some important points. In the geometric method [HI IE] , the 
authors start with a four dimensional harmonic oscillator V4 = | J2t=i u2n a s a an d then 
they obtain the nonlinear system by reducing the dimension from n = 4 to n = 2. In 
the method presented in this paper we also start with an harmonic oscillator but then we 
obtain the nonlinear system by a deformation of the complex quadratic equation. Thus, in 
both cases the starting point is the (linear) harmonic oscillator, which is superintegrable 
in any dimension. 

Finally, it is natural to think that the introduction of the deformation (functions hi(x{)) 
in the quadratic equation can also be applied to equations of order higher than 2. There- 
fore, a natural generalization of this formalism would be the search of new superintegrable 
sytems by introducing deformations in higher powers K™ with m > 2 of the functions K{. 
Note that we have restricted the study to deformations of the real part of the functions 
Kf; nevertheless in the more general case the deformations could be introduced also in 
the imaginary part of the complex functions. This question, as well as some other related 
problems (as the properties of the quantum version of this system) are open questions to 
be studied. 
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